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Abstract: It is shown that in some multi-supergraviton models, the contributions to the effective 
potential due to a non-trivial topology can be positive, giving rise in this way to a positive 
cosmological constant, as demanded by cosmological observations. 

hj ; 

Os '• 

^ ! 1 Introduction 

(N '■ 

Renewed interest in the study of multi-graviton theories [Q owes, in particular, to the fact that 
these formulations resemble higher-dimensional gravities in the presence of discrete dimensions. 
These classes of discretized Kaluza-Klein theories are now in fact under the focus of attention 
due to their primary importance for the realization of the dimensional deconstruction program 
[3] . Moreover, multigravitons can be also related with discretized brane- world models [I]. 

In spite of the absence of a consistent interaction among the gravitons, one can think of 
possible couplings in the theory space. In particular, in a recent paper [Sj, a multi-graviton 
theory with nearest-neighbor couplings in the theory space has been proposed. As a result, a 
discrete mass spectrum appears. The theory seems to be equivalent to Kaluza-Klein gravity 
with a discretized dimension. 

In a previous paper concerning multi-graviton theory |2J, we have shown by means of an 
explicit example, namely a discretized Randall-Sundrum (RS) brane-world jjj, that the induced 
^ ■ cosmological constant becomes positive provided the number of massive gravitons is sufficiently 

large. 

In the present paper, we would like to show that an alternative mechanisms can also give rise 
to positive contributions to the cosmological constant. In particular we shall consider a multi- 
supergraviton example with few gravitons, in a manifold (bulk) with non trivial topology. We 
shall show that in such a model a positive cosmological constant A can be generated, due to the 
presence of positive topological contributions. Moreover, by a suitable choice of the topological 
parameters, the number obtained for A can reach a value perfectly in accordance with result 
obtained from recent cosmological observations [Sj. 
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2 The multi-graviton and multi-supergraviton models 
2.1 The graviton model 

We start by considering the Lagrangian for the spin-two field h^ v with mass m 
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C m = C -^-(h^ -h 2 ) -2(mA^ + d^)(d"h^-d fl h) 
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1 
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(0 M A„ - dvAJ (d^A v - d v A») , (2.1) 



where Cq is the Lagrangian of the massless spin-two field (graviton) h^ u (h = ' ) 



Co = -\d x h^d x h^ + d x h\d u h^ - d^d u h + l -d x hd x h , (2.2) 



while A^ and (p are Stiickelberg fields [Sj. 

The multi-graviton model is defined by taking iV— copies of (j'2.1j) with graviton h nfll , and 
Stiickelberg fields A n ^ and y? n . Here, we propose a theory defined by a Lagrangian which is 
taken to be a generalization of the one in 0. It reads 



* = £ 



N-l r j 



n=0 



■zdxK^frhg + d x h*d u K v - djtfdvK + -d x h n d A h n 



~ (m 2 Ah npLV Ah^ - (Ah n ) 2 ) - 2 (mAU£ + d*V n ) (5"/^ - d„h r 
~(d lt A w -d v A nil )(& t AZ-ff'A!fi . 
The A and A 1 " are difference operators, which operate on the indices n as 



(2.3) 



N-l N-l N-l 

A(f) n = ^2 a k<Pn+k , AT0 n = Ofc^n-fc , J! = ' ( 2 ' 4 ) 

fc=0 fc=0 k=0 

where the are N constants and the N variables 4> n can be identified with periodic fields on 
a lattice with N sites if the periodic boundary conditions 4>n+ N = 4>n are imposed. The latter 
condition in (|2.4JI assures that A becomes the usual differentiation operator in a properly defined 
continuum limit. 

The eigenvalues and eigenvectors for A are given by 4 

N-l 

A<% = i^ft , Vip=£ a n e 2mn ^ N , (2.5) 



n=0 

g2irinp/N 

„ =-. p = 0,l,2,3,... (2.6) 



'N 

By using ()2.4|) in the latter equation and assuming a n to be real one gets the relations 

^o = 0, Hp = —[iN-p , HN-p = , (2-7) 



4 Please note that here we use a different notation with respect to one used in Refs. [§] and In fact, in 

order to avoid confusion with masses, we have replaced the eigenvalue m with /i and the index M with p. 
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which, for any fixed N, permit to obtain the whole spectrum of the theory. 

Then we see that the Lagrangian l|2.«-ij) describes a massless graviton and N — 1 massive 
gravitons, with masses M p = \fi p \ (p = 1, 2, . . . , N — 1). It must be pointed out that the massive 
gravitons always appear in pairs which share a common mass and, moreover, the complex mass 
parameter \i v can be arbitrarily chosen, just by properly selecting the coefficients in 1)2,51) 6 . 
As discussed in [3], the multigraviton model can be regarded as corresponding to a Kaluza-Klein 
theory where the extra dimension lives in a lattice. 

As a specific example, we now consider the two-brane Randall-Sundrum (RS) model (for 
a recent review see [1U] ^) . In this model, the masses of the Kaluza-Klein modes are given by 

M p = ^ , z c = l (W z - l) , p = 0, 1, 2, ... (2.8) 

where I is the length parameter of the five-dimensional AdS space and 7rr c the geodesic distance 
between the two branes. 

Motivated by this last equation (|2,8|) . we consider an A = 2N' + 1 graviton model, with the 
graviton masses being given by 

[ f> p = 0,l,..-,iV', 

M P S _ I {N-p 1 p = N > + hN > + 2i ...,N-l = 2N'. 

Those are solutions of Eq. ([2.5(1 . with the choice ao = and, for any n > 1, 

• 2tt N'n\ ■ 2tvu 



Im 



\ — g 2N' + 1 ] g 2JV' + 1 



(2JV' + l)z, 



g 2JV' + 1 



(-1)» 27T ^ 



A^z c sin (^) 



(2.9) 



We see that plays here the role of a cutoff of the Kaluza-Klein modes. 

In previous models of deconstruction IE] , mainly nearest neighbor couplings between the 
sites of the lattice have been considered. As a consequence, on imposing a periodic boundary 
condition, the lattice then looks as a circle. Departing from this standard situation, in the model 
considered here we have introduced non-nearest-neighbor couplings among the sites. That is, a 
site links to a number of other ones in a rather complicated way. In this sense, the lattice in 
the present model is no more a simple circle but it looks more like, say, a mesh or a net. Let us 
assume that the sites on the lattice would correspond to points in a brane. If the codimension of 
the spacetime is one, the brane should be ordered, resembling the sheets of a book. One brane 
can only interact (directly) with the two neighboring branes. However, if the spacetime is more 
complicated and/or the codimension is two or more, the brane can directly interact with three 
or more branes, an interaction that will be perfectly described by our model corresponding to 
this case. For example, a site on a tetrahedron connects directly with three neighboring sites. In 
this way, the non-nearest-neighbor couplings we here consider may quite adequately reflect the 
structure of the extra dimension. In this respect our model is very general and opens a number 
of interesting possibilities. 
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2.2 The supergravity case 



By using the same sort of techniques described above, the multi-graviton model can be gen- 
eralized to the supergravity case, just by starting with a supergravity theory in 5-dimensions 
and implementing deconstruction by way of replacing the fifth spacelike dimension with a one- 
dimensional lattice containing iV-points. A multi-supergravity model of this kind has been 
proposed in Ref. to which the interested reader is addressed for details. Here we shall only 
write down the essential aspects which will be used in what follows. 

In the 5-dimensional linearized supergravity theory, the number of bosonic degrees of freedom 
is 8, 5 due to the massless graviton and 3 due to the massless vector (gauge) field and the number 
of fermionic degrees of freedom is 8 too, due to the the complex Rarita-Schwinger field (4x2). 

The deconstruction process now consists in replacing the fifth dimension in the action of 
spin two+vector+Rarita-Schwinger fields with N— points and the derivatives with respect to the 
corresponding variable with the operator A as in Eq. (|2,5|) . In this way one gets a complicated 
action in 4 dimensions, similar to the one in Eq. (|2.3j) . but with vector and fermion parts too. 
It contains a spin-2 field hav (the graviton), but also scalar, vector and fermionic fields. More 
precisely, in the massless sector one has 8 degrees of freedom due to bosons (graviton (2 d.o.f.), 
gauge and Stiickelberg vectors (2+2 d.o.f.), a Stuckelberg scalar and the fifth component of the 
gauge field (1+1 d.o.f.) and 8 degrees of freedom due to fermions (complex Dirac and Rarita- 
Schwinger fields), while in the massive sector one has again 8+8 degrees of freedom, but only due 
to a massive graviton, vector and Rarita-Schwinger fields. As in the pure-gravity case, one has 
N copies of such fields and their masses — obtained by imposing periodic boundary conditions — 
are always given by means of Eq. ()2.5|) , that is 

N-l 

A< = i^ft , On e 2mnp/N , (2.10) 

n=0 

„2irinp/N 

€ = €+N = nx-- P = 0,1,2,3,... (2.11) 



N 

On the other hand, for fermion fields anti-periodic boundary conditions could also be considered. 
In such case one gets a different spectrum, given by means of the following equations 

N-l 

= ifaft , ifip = E °n e^CP+V*)/" , (2.12) 



n=0 

e 2mn{p+l/2)/N 



— ~Vn+N — ' p — 0,1,2,3, ... (2.13) 

It has to be noted that with boundary conditions of this sort, there are no massless fermions 
and this is a consequence of the explicitly breakdown of global supersymmetry. 



3 The induced cosmological constant 

We now turn to the evaluation of the induced cosmological constant for the N— graviton and 
super-graviton models discussed in the previous section. To this aim — the main one in the 
present paper — we shall compute the one-loop effective potential using zeta-function regular- 
ization |12 | I13|: needless to say, other regularization schemes could work as well. First of all, we 
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compute the effective potential for a free scalar field with mass M, since this corresponds to the 
contribution of each degree of freedom to the one-loop effective potential of our theories. 

In the zeta-function regularization method, the one-loop contribution to the effective poten- 
tial is given by 

v e}} = -^c'(oiiv^) = ~<am - -^c(oiL)io g ^ 2 , (3.i) 

V being the volume of the manifold and C( s l-^) t ne zeta function corresponding to the Laplacian- 
like operator L = — A 2 + M 2 , with M a positive constant. The arbitrary parameter /i has to be 
introduced for dimensional reasons. It will be fixed by renormalization at the end of the process. 

The manifold we are considering in the present paper is a flat one with non trivial topology 
of the kind M = Ex T 3 . The simplest case M = M 4 has been already considered in [51 ITT]. 

The operator L has the form 

L = --^+L 3 , L 3 = -A 3 + M 2 , (3.2) 

A 3 being the Laplace operator on T 3 . The zeta-function is expressed in terms of the heat trace 
via the Mellin representation. The heat traces read 

Tie- tL = VJC(t\L), Tre- ti3 =V 3 JC(t\L 3 ), K(t\L) = K( - t ^\ (3.3) 



'4mt 



As a result 



1 roo \r poo 

«.|l) = W) j o dtt -i Tre -* = 7 =-l At .-v* mL3) 

- V ^-^«s-m\L 3 ), (3.4) 



47rr(s) 

<^(s — I/2IL3) being the zeta-function density on T 3 and V3 = (2-7rr) 3 the "volume" of the torus 
with "radius" r. The heat kernel and zeta function on T 3 are well known. In the Appendix 
lAl for the reader's convenience, we summarize some useful representations that will be used in 
what follows (for a review, see |LS|). 

Using expressions ()3.4|) and (jA. 1|) one realizes that the zeta function can be written as the 
sum of two terms, that is 

C(s\L) = ftHL) + (t(s\L) , (3.5) 

where £0 is the same one has on H 4 , namely 

VT(s - 2)M 4 - 2s VM*- 2s 
C ° (S|L) " 167r 2 r( S ) -l6vr 2 ( S -l)( S -2)' (3 - 6) 

while represents the contribution due to the non-trivial topology, which explicitly depends 
on the topological parameter r. Expression ()3.6|) is also the leading contribution to the whole 
zeta function in a power series expansion for large values of M. 
Recalling now (|A,4|) . we obtain 



, , ,^ VT(s — 3/2) cos tts Af 4_2s , w 2 iN o /2 , 

(t(s\L) = i ' ' - / du G(Mru) (u 2 - lf' 2 ~ s . (3.7) 



8vr 5 /2r( s ) A 
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Observe that the topological contribution vanishes at s = 0, and this means that 

C(0|L) = Co(0|X) = . (3.8) 

Using Q3.1J) . for the one-loop effective potential we finally have 

It is interesting to note that for scalar fields, in the large mass case the topological contribution 
is always negative, and it is negligible with respect to the standard Coleman- Weinberg term. 

As we have anticipated above, the parameter \x has to be fixed by a renormalization condition. 
To this aim, here we follow Ref . ^1] . The total one- loop effective potential is of the form 

V eff = V R (i,) + V e f f (v), (3.10) 

Vr(/x) being the renormalized vacuum energy. For physical reasons, the last expression has to 
be independent of /x, and this means that 

,^o, (,n, 

from which it follows that 

y^) = ^(^) + — log^f, (3.12) 

being the renormalization point which has to be fixed by the condition V r {^jlr) = 0. In this 
way, we finally get 

M 4 / M 2 3 \ 

M 4 f°° , „. w o N o /9 M 2 ^ K 2 {2nMr\n\) 

ne'E ;n^0 

Vr(r) represents the contribution coming from the non-trivial topology, which for scalar fields 
is always negative. We also note that, as a function of the topological parameter r, Vr(r) can 
reach, in principle, any negative value. In Eq. (|3.14j) . K v are the MacDonald's (or modified 
Bessel) functions. 

Before proceeding with the computation of the induced cosmological constant corresponding 
to the models we have discussed in Sect. 2, we first analyse here the behavior of Vr(r) as a 
function of r. To this aim, we consider the two different regimes Mr <C 1 and Mr 1. 

For the case Mr <C 1, using ()A.5|) in 1)3. 14j) we get 

V T (r) = -tt^ / duG{Mru){u 2 - if' 2 

oo 

dxG{x)(x 2 -MV) 3/2 

Mr 



12vr 2 Ji 
1 

12vr 2 ? 
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12vr 2 r 4 



dx 1 



Mr 



M 2 r 2 



3/2 



/*1 /*oo 

+ y dxG (x)x 3 + J dxG(x)x 3 + 0{M 2 r 2 ) 



1 



32vr 2 r 



2».4 



E 

73 



1 



+ 0(MV) 



neS ;|n|^0 



1 



64vr 2 r 4 



+ 0(M 2 /r 2 



(3.15) 



We thus see that in this limit the leading term does not depend on M, and that it can be 
arbitrarily large, with a suitable choice of the parameter r. The series in the latter equation has 
been computed numerically. 

On the contrary, in the opposite regime, Mr 3> 1, using (|3.14|) and the asymptotic expansion 
for the Bessel function, we obtain 

M 2 

V T {r) = 



16vr 4 



E 

3M 4 

e 



K 2 (2irMr\n\ 



2nMr 



32^ 4 (Mr)5/2 - +••• ( 3 - 16 ) 

In this limit the topological contribution could indeed be arbitrarily small. In Fig. 1 the whole 
behavior of the topological part of the effective potential is drawn. In order to work with 
dimensionless variables we have introduced the function Vr(y) = r 4 Vr(r) of the dimensionless 
variable y = Mr. The graphic corresponds to the exact expression for Vr(y), as given e.g. by 
the first lines of Eq. (|3.15^ , multiplied by 3 • 6Air 2 . A very smooth transition from the behavior 
corresponding to Mr <C 1, Eq. (|3.15() . to the one for Mr 3> 1, Eq. (|3.16|) . is revealed. In Fig. 2, 
the corresponding graphic of the full effective potential, Eq. (|3.13|) . is drawn, again as a function 
of y and setting [i^r = 1. 

At this point, the effective potential —and, as a consequence, the induced cosmological 
constant for the models we are interested in— can be obtained by adding up several contributions 
of the kind (ETHl . 



3.1 The multi-graviton model 

We start with the explicit example of multi-graviton model given by (|2.9|) . in which there is a 
single massless graviton and (N — l)/2 couples of massive gravitons, with masses given by 

7TO iV — 1 

M o = |/i O |=0, M p = \fi p \ = -^, p = 1,2,...,—— . (3.17) 

z c z 

On the manifold M = 1R x T 3 , the massless graviton gives no contribution to the effective 
potential, while it does appear explicitly on manifolds with a non-vanishing curvature. Since 
the massive gravitons always show up in pairs, in order to perform the computation of the 
effective potential, it is sufficient to consider only one half of the whole massive spectrum. 
Moreover, we have to take into account that each massive graviton contributes with five scalar 
degrees of freedom. After these considerations have been properly taken into account, for the 
effective potential we get the following expression 

■V-l 

s M 4 / M 2 3\ 
V eff = 10 V — ln-f-- 
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V T (y) 




One can see that, as for the non-compact flat case (see Ref. [Sj for details), in order to have a 
(small) positive induced cosmological constant one has to consider a large value of N, that is, 
a huge number of massive gravitons. In this respect, the torus topology does not improve the 
situation. As we are now going to show, this is no longer the case for the multi-supergraviton 
model. 



3.2 The multi-supergraviton model 

Here we have to distinguish two cases: the first one corresponds to the choice of periodic bound- 
ary conditions in both the bosonic and fermionic sectors. In such situation, the degrees of 
freedom due to bosons exactly compensate the degrees of freedom due to fermions. Moreover, 
for any massive boson there is a fermion with the same mass and, since the contribution to the 
effective potential of any fermionic degree of freedom is opposite to the contribution of a bosonic 
degree of freedom, it turns out that the induced cosmological constant vanishes, independently 
of the mass spectrum. 

In the second case, that is when anti-periodic boundary conditions are imposed in the 
fermionic sector, the situation changes completely, since the fermionic mass spectrum becomes 
really different with respect to the bosonic one. For example, by choosing N = 3 the 
solutions of Eqs. 1)2. lOj) and lj2.12j) give 

M = , Mi = M 2 = m , for bosons , (3.19) 
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2.5' 



Figure 2: The exact expression for V e ff{y) = r4 K//(j), Eq. Q3.13|) . as a function of y = Mr. 



M = M 2 



rn 



V3' 



Mi 



2 m 
73 



for fermions , 



(3.20) 



m being an arbitrary mass parameter. 

By taking into account the number of degrees of freedom, the one-loop effective potential 
becomes, in this case 



V 



eff 



4^2 



In 



AM? 



M 4 



4tt 2 y 

8vr 2 
,16 



, M 3 \ 



Mf( ln W_l 



du G{M ir u) (u 2 - 1) 3/2 

du G{M ru) (u 2 - 1) 3/2 

du G(M iru ) (u 2 - l) 3 / 2 



4M 4 




3tt 2 


Jl 


) + 


2Mf 
3vr 2 



, 2 1 
36^ l0g 3^ + yT 



(3.21) 



Vt being the sum of all the topological contributions. As one sees, the first term on the right- 
hand side of the latter equation is always negative, but the whole effective potential could be 
positive due to the presence of the topological term. For example, in the regime mr -C 1, from 
(|3.15j) one has 



V T 



1 



Veff > 



(l 2 lfr 
for mr < - log — 77- 
V 9 3 9 



-1/4 



1.4, 



(3.22) 
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while in the opposite regime, mr 3> 1, by using (|3.16j) one can see that the topological contri- 
bution although still positive it is negligible, and thus the effective potential remains negative. 

In Fig. 3, the corresponding graphic of the full effective potential, Eq. (|3.21|) . is drawn, 
again as a function of y = mr. 



V eff (y) 




0.8 1 1.21.4 



Figure 3: The exact expression for V e ff(y) = r 4 V e //(r), Eq. (|3.21() . as a function of y 



: mr. 



4 Conclusions 

In this paper, we have computed the effective potential corresponding to a multi-graviton model 
with supersymmetry in the case where the bulk is a flat manifold with the topology of a torus 
(more precisely Ex T 3 ), and we have shown that the induced cosmological constant could be 
positive due to topological contributions. 

In previous papers [SI multi-graviton and multi-supergraviton models have been con- 
sidered in IR 4 . It has been shown that in the multi-graviton model the induced cosmological 
constant can be positive, but only if the number of massive gravitons is sufficiently large, while in 
the supersymmetric case the cosmological constant can be positive if one imposes anti-periodic 
boundary conditions in the fermionic sector. Note that the topological effects discussed above 
may also be relevant in the study of electroweak symmetry breaking in models with a similar 
type of non-nearest-neighbour couplings, for the deconstruction issue |15j . 

In the case of the torus topology, the topological contributions to the effective potential have 
always a fixed sign, depending on the boundary conditions one imposes. In fact, they are negative 
for periodic fields and positive for anti-periodic fields. This means that the torus topology 
provides a mechanism which, in a most natural way, permits to have a positive cosmological 
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constant in the multi-supergravity model with anti-periodic fermions, being the value of such 
constant regulated by the corresponding size of the torus. 1 In this situation one can most 
naturally use the minimum number, N = 3, of copies of bosons and fermions. 

We finish with the remark that it would be interesting to apply the deconstruction scheme of 
Ref. |S] also for the case of two latticized extra dimensions, which in the continuous limit would 
contain the orbifold singularity. This analysis might have a quite interesting impact on brane 
running coupling calculations |17j . 
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A Zeta function on the torus 

Eigenvalues of the Laplacian on the 3-dimensional torus are of the form A n = n 2 , n £ 'S? , and 
thus the corresponding heat kernel is given by 

p -tM 2 p-tM 2 

K(t\L 3 ) = G — £ e-^ 2 = £ e -W/« , (A.l) 

being V 3 = (2-7rr) 3 the "volume" of T 3 . Using the expression above, the zeta function of this 
Laplacian can be written as 

((s\L 3 ) = ( (s\L 3 ) + ( T (s\L 3 ) , (A.2) 

where the contribution 

t{s ir ^ _ - 3/2)M 3 / 2 - 2 * 

C ° (S|L3) " (4vr) 3 / 2 I» ' (A - 3) 

comes from the n = term and it is the same one has for 1R 3 , while Ct corresponds to the 
contribution due to the non-trivial topology. Such term can be written in different ways, for 
instance, as an infinite sum of Bessel functions. 

In Refs. |12j one can find many interesting results concerning zeta functions and heat kernels 
corresponding to operators on manifolds with constant curvature. In particular, on the torus 
one has the very nice representation 



Ct(z\L 3 ) = M ^ Z {l S ^ z) Z ] x du G(Mru) (n 2 - 1)1- 
n z - 2 ^ /Af\ 3 / 2 - 2 



4T(z) ^ \r\n\ 



£ b^T K 3/2 _ z (27rMr\n\) , (A.4) 



X A more crude analysis for the pure scalar case already hinted towards this conclusion. However, the sign issue 
was there not easy to fix (16) . the reason being now clear. 
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where G(x) is given by 

G( x )= Yl e ~ MAx = + 4 E ' 



7T 2 ^3 (n 2 + X 2 ) 2 



-1 + -^t + ^ G 0( X ) , ( A - 5 ) 



7T 2 X 3 7T 2 



Gq(x) being the regular function 
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